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ABSTRACT 

We study the star formation history for a sample of 154 galaxies with stellar mass 
10 10 < M* < 10 12 Mq in the redshift range 0.7 < 2 < 0.9. We do this using stellar popu¬ 
lation models combined with full-spectrum fitting of good quality spectra and high resolution 
photometry. For a subset of 68 galaxies (M* > 10 11 Mq) we additionally construct dynamical 
models. These use an axisymmetric solution to the Jeans equations, which allows for velocity 
anisotropy, and adopts results from abundance matching techniques to account for the dark 
matter content. We find that: (i) The trends in star formation history observed in the local 
universe are already in place by z ~ 1: the most massive galaxies are already passive, while 
lower mass ones have a more extended star formation histories, and the lowest mass galax¬ 
ies are actively forming stars; (ii) we place an upper limit of a factor 1.5 to the size growth 
of the massive galaxy population; (iii) we present strong evidence for low dark matter frac¬ 
tions within 1 R e (median of 9 per cent and 90th percentile of 21 per cent) for galaxies with 
M* > 10 11 Mq at these redshifts; and (iv) we confirm that these galaxies have, on average, a 
Salpeter normalisation of the stellar initial mass function. 
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1 INTRODUCTION 

Currently, one of the key quests of Astrophysics is to under¬ 
stand and model the processes that guide the formation and evo¬ 
lution of galaxies. Great strides have been made over the past 
few decades and with the advancement of technology, such as 
ever larger telescopes taking ever larger surveys of thousands 
of ga l axies within an ever larger redshift range ( e.g. 


these quiescent galaxies have undergone a size evolution; with the 


200C; JCollessrtalJJ^OOlJ^ LeFevreetalJ_ 20051 : iNewman et al 
2013 ; ICoil et all 2011 ; Kochanek et alJ 2012h. the advent of new 


York et al 


techniques such as g ravitational lensing le.g. iBolton et alJl2006l : 
iHevmans et alJ l2012h. and galaxy surveys using integral field 
spectroscopy (e.g, IdeZeeuwgtalJ 2002|^ Caggehan^dJ 1201 lal : 


ISdnchez et alJl2012i:IScott etal 20irTBundvetaFl 201 5t) . 

Recent observational evidence suggests that the star forma¬ 
tion rate of the universe peaked at z ~ 2.5 and that by 2 ~ 
1 half of the stellar mass of the universe today was already 
in place i Bundv. Ellis & Conselicel 120051 : iMortlock et alJ 1201 ll : 
iMuzzin et al. 2013h . The decreasing star formation rate, referred 
to as quenching, is mass dependent with the more massive galax¬ 
ies being quenched earlier. Also, the comparison of the most mas¬ 
sive galaxies (A/* > 10 11 ATq) at high and low redshifts show that 
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dated with minor mass growth, suggesting that these growths may 
be driven by minor merger where the size of the galaxy grows to 
the second power of the added mass through virial arguments, un¬ 
like maj or mergers where the size grows l i nearly to the increase 
injnassjNaab Ji Johansson_&_Ostriker 2009; Bezanso n et al] [20091 : 
van der Wei et alj |2009i) . Additionally, recent works have pointed 
out that a significant part of the observed size growths in the popula¬ 
tions of quiescent galaxies, especially at lower masses, may be due 
to progenitor bias, wherein the addition of large recently quenched 
galaxies contrib ute to the observed increase in the mean s ize of the 
population fe.g. lCarollo et al.ll2013l : IPoggianti et al.ll2013h . 

Regardless of what the process for the growth of the galaxy 
size, and its stellar mass may be, there is strong evidence in¬ 
dicating that, for the most massive galaxies, most of the ad¬ 
ditional stellar mass is added to the outskirts of the galaxies, 
while the central regions remain mostly unperturb ed (Trujillo et ah 
120061 : IBezanson et al] 120091 ; Ivan Dokkuni et all l20ld) . The end 
result of this merging process are the most massive galaxies 
in the nearby Universe which are found to be slowly rotating 
l lEmsellem et al.ll201ll ; ICappellari et ~ai] |2013al) . they have cores in 
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their surface brightness pro files (Faber et all 1 19971 : iFerrarese et alJ 
l2006l :lKramovic et al.||201 3), and are embedded in extended stellar 
envelopes (Kormendy et al. 200^). 


The situation appears radically different for less massive 
(M*<10 11 Mq) passive galaxies. At the present day, they are 
structurally differen t, and appear to hav e followed a differ¬ 
ent evolution path ( 


metric (iKrajnovic et al. 

2011), they contain disks like spiral 

galaxies (Cappellari et al 

201 lb) and are dominated by rotation 

(Emsellem et al. 201 lj; Cappellari et al. 2013a). These fast rotating 


galaxies follow the same mass-size relation, and have the sa me mas s 
distribution, both in dense clusters as in the field ICappellarfcO 13l) . 
indicating they experienced an insignificant amount of merging dur- 
ing their evolution, in a greement with redshift evolution studies 
Ivan Dokkum et allhold) . 

Due to the recent advances in the techniques of stellar popula¬ 
tion modelling and redshift surveys, a key addition to this emerging 
picture of galaxy evolution is provided by studies of the stellar pop¬ 
ulations of galaxies throug h cos mic time. 

The work of Schiavon_et aljj 20061) , using spectra from the 
DEEP2 sur vey i Newman et alJ 2013l) ~compared to local SDSS 
lYork et al.M200C 1 results, suggests that the evolution of the red- 
sequence galaxy population is not consistent with a passive evo¬ 
lutionary model. Instead, they propose that the red-sequence pop¬ 
ulation should either continue to host some level of star formation 
(“Frosting”) to present day or have newly quenched galaxies joining 
the red-sequence galaxies between 2 ~ 0.9 and today. 

IChoi et al] ( 120141) study quiescent high redshift galaxies via a 
full spectrum fitting of stacked galaxy spectra to derive the stellar 
ages and elemental abundances of Fe, Mg, C, N and Ca. The work 
uses optical spectra of local gal axies taken from the SDSS and spec¬ 
tra from the AGES dKochanek et al.ll2012l) survey within a redshift 
range of 0.1 < 2 < 0.7. They find negligible evolution in elemen¬ 
tal abundances at fixed stellar mass. For the most massive galaxies 
they measure an increase in stellar age consistent with passive evo¬ 
lution since 2 0.7. While at masses below ~ 10 10 ' 5 Mg, the 

data permit the add ition of newly quenched galaxies. 

iGallazzi et ali (120141 Hereafter G14) study a sample of 70 qui¬ 
escent and star-forming galaxies at 2 tn 0.7, above a stellar mass 
of 10 10 Mg. They derive the stellar age-mass relation of the galax¬ 
ies, which they compare with the one derived in a similar manner 
in the local universe. They find that taken as a whole, passive evo¬ 
lution cannot represent the evolution of galaxies in the last ~ 6 
Gyr. In fact, although the shape of the stellar age-mass relationship 
between the two redshifts is similar, the offset is inconsistent with 
passive evolution. This is agreement with their observed metallic- 
ity differences with redshift. They propose a mass-dependent star 
formation history (SFH) to explain the observations. 

Here we use full-spectrum fitting to explicitly determine trends 
in the star formation history of a sample of 154 galaxies at 0.7 < 
2 < 0.9. Furthermore, we investigate the correlation between the 
stellar population and the physical parameters of the galaxies. We 
also present results on the dynamical modelling of a subset of 68 
galaxi es. This subsample is the same we analysed in our previous 
work Ishettv & CannellarillO14) . where we studied the Initial Mass 
Functi on (IMF) mass normalisation and concluded it is consistent 
with a lSalnete'i] (1955) slope. Here, we improve upon the dynam¬ 
ical models by accounting for the dark matter of the galaxies via 
abundance matching techniques. 

In Section 2 of the paper, we describe the observational data 
that we use within this study while in Section 3 we discuss the se¬ 
lection criteria that have been implemented in the course of this 


analysis, along with a comparison of our galaxy sample with the 
parent sample. In Section 4, we present the various methods used to 
analyse the data followed by a discussion of the results in Section 5. 
In Section 6, we provide a summary of the results. 

In this study, we assume a flat universe with the following cos¬ 
mological parameters: D m = 0.3, Da = 0.7, Ho = 70 km s -1 
Mpc -1 . 

2 DATA 

2.1 Spectral Data 

The spectral data used in this study was taken from the DEEP2 sur¬ 
vey. The survey is an apparent magnitude limited, Rab = 24.1, 
spectroscopic redshift survey of 4 fields covering an area of 2.8 
deg 2 across the sky characterised by low extinction. The observa¬ 
tion targets for the survey were selected using an algorithm based 
on the Canada-France-Hawaii Telescope ( CFHT) obs erved object 
size, luminosity, (B-R) and (R-I) colours ( iCoil et all 120041) and a 
surfac e brig htne ss cu t-off which eliminated wrongly classified ob¬ 
jects iNewman et alJl2013h . Of the four fields o f DEEP2, only the 
Extended Groth Strip (EGS: [Groth et alJ (1 1994) ) has the required 
high resolution HST data to construct dynamical models and accu¬ 
rately measure galaxy siz. Hence this study is restricted to this field 
of the DEEP2. 

The observations for the survey were done using the DEIMOS 
spectrograph on the Keck-2 telescope, covering an observed wave¬ 
length range of 6,100A - 9,100A. The highest resolution grating, 
1,200 line mm' , was used to achieve a spectral resolution of R 
~ 6, 000 within the l"x ~ 7" slitlets of the instrument. The ob¬ 
servations were taken in 3 exposures of 20 minutes each, with the 
FWHM generally ranging between 0" .5 and l".2 due to seeing. 

The reduction of the 2D spectrum of each galaxy and the ex¬ 
traction of its ID spectra was done by the DEIMOS spec2d pipeline. 
In our study we use the ID spectra extracted via the Boxcar tech¬ 
nique from the central ~ 1" x 1" window. Due to the design of 
the survey, the observed ID spectra is divided into 2 halves along 
wavelength. Here we use the blue half of the two, since the CaK 
(3,934A) and G band (4,304A) stellar lines fall within this half for 
galaxies in the redshift range that we are studying. 


2.2 Photometric Data 

Our photometric data was taken from the AEGIS data product 
1 Davis et al.ll200 7). specifically from the Hubble Space Telescope 
(HST) GO program 10134 (P.I.: M. Davis). These are HST /ACS 
images taken in the F606W (V band) and F814W (I band) filters. 
The observations are limited to 5 o measurements of Vf 606 W = 
26.23 (AB) and 7fsi4W = 25.61(7173) within a circular aperture 
of 0”3. This program covered 0.197 deg 2 of the EGS, overlapping 
entirely with the region observed by DEEP2. The observations were 
made in a four point dither pattern, which was processed by the 
STSDAS MultiDrizzle package to produce a final mosaic of the 
field with a pixel scale of O'.' 03. In this work, we use the F814W 
photometry as this translates to the vega B band at the redshifts that 
we are studying. 


3 GALAXY SAMPLES 

In this study we use two samples extracted from the large DEEP2 
survey catalog. The first one consists of 154 galaxies and was used 
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Galaxy Sample vs Parent Sample 



Figure 1 . In this plot, we compare our galaxy samples: 68 galaxies in our secondary sample (red points) and 154 galaxies in our full sample (blue and red 
points), against the parent sample of 2, 896 galaxies in the same redshift range and field from the DEEP2 survey (orange points). In the left panel, we present 
the colour-magnitude plot where we can distinguish between the red sequence and blue cloud of the parent sample. It is evident that our galaxy samples have a 
measurably higher magnitude cut-off, compared to the parent sample, which is due to our selection criteria. Within this higher magnitude cut-off, our full galaxy 
sample provides a good representation of both the red sequence and the blue cloud. We also notice that our secondary galaxy sample is a good representation 
of the red sequence, while also including a few very luminous galaxies of the blue cloud. If we assume that the Rab magnitude is a tracer of stellar mass, the 
right panel can be approximated as a stellar mass-size plot for the galaxies. We see here that our galaxy samples mainly consist of galaxies above a certain 
surface brightness. 


to measure trends in the galaxies star formation histories via full- 
spectrum fitting. The second one is a subsample of 68 galaxies of 
the first, for which we construct dynamical models. Details of the 
adopted selection criteria are presented below. 


3.1 Selection Criteria 

Our first selection criterion is based on the quality flags provided 
by the DEEP2 survey catalogue, i.e. we selected galaxies that have 
either “secure” or “very secure” redshift estimation based on the au¬ 
tomated and visual spectral redshift estimation by the DEEP2 team. 
We select galaxies with redshift 0.7 < z < 0.9. In this study we 
have only used galaxies in the EGS field due to the lack of availabil¬ 
ity of high resolution HST photometry in the remaining fields. By 
extracting 9" x 9” thumbnails, from the AEGIS mosaics, centred 
at the galaxy coordinates provided by the DEEP2 catalogue, we vi¬ 
sually confirm the presence of photometry of the galaxies. These 
thumbnails are later used to parametrize the photometry of the 
galaxies and further details on this is provided in Section l4~3l From 
a total of 43,569 galaxies in the DEEP2 Galaxy Redshift survey 
1,249 galaxies passed this stage of the selection process. Of these, 
1,240 (99.3%) galaxies have an apparent Jf 814 w < 24(AB). 

The second selection criterion is based on the quality of the 
galaxy spectra. We calculate the mean signal-to-noise ( S/N ) of 
a galaxy spectrum as the median flux in the galaxy spectrum di¬ 
vided by the standard deviation of the residuals of the spectral fit to 
that spectrum. In this instance, the spectral fitting was done with a 
set of empirical stellar spectral library, while masking gas emission 
features. Further details on this spectral fitting is presented in Sec- 
tion l4.ll We only select galaxies with S/N > 3. We further elimi¬ 
nate galaxies that do not appear to have significant absorption fea¬ 
tures, i.e. clear spectral absorptions which stood out visually above 


the noise. This selection criteria reduced the galaxy sample to 175 
galaxies. 

Our third selection criteria is based on the visual inspection 
of the thumbnails of the galaxies, taken from the F814W images 
by the HST/ ACS. We have removed 21 irregular galaxies from our 
final sample as these galaxies are likely to be disturbed and would 
confuse the trends being investigated in this study. The selection cri¬ 
teria based on the galaxy spectra and imaging quality, dramatically 
reduces our galaxy sample to 154 galaxies. 

In this study, we also create dynamical models for our galaxy 
sample, however this requires some additional selection criteria. To 
create the dynamical models we need to parametrize the photome¬ 
try of our galaxies to high accuracy within the central regions from 
which the stellar kinematics are observed. We exclude galaxies that 
have disturbed photometry, such as dust lanes, or non-axisymmetric 
features, and hence this galaxy subsample is biased towards early- 
type galaxies with smooth light profiles. Further details on this se¬ 
lection criterion is given in Section POI 

Also, the dynamical modelling of galaxies assume a spatially 
constant stellar mass-to-light ratio ( M/L ), however this assump¬ 
tion becomes inaccurate for galaxies with multiple significant star 
formation events. Therefore, after deriving the SFH of our galax¬ 
ies, we eliminate galaxies that require more than one significant 
star formation episode to reproduce their spectrum. In addition to 
this, the ATLAS 3D group found that young galaxies, as inferred 
by their strong H/3 absorption feature, tend to show strong gradi¬ 
ents in the stellar M/L. For t his reason, these galaxies whe re re¬ 
moved from their IMF studies (Cannellari et alJEoTMlTQPbl . here¬ 
after C 13b). In this study, we use full-spectrum fitting to the galaxy 
spectra to identify the best fitting single stellar population (assum¬ 
ing solar metallicity) of the galaxies. We then classify any galaxy 
with the best fitting population of 1.2 Gyrs or younger as a young 
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galaxy and eliminate them from our secondary sample. This nar¬ 
rows our secondary galaxy sample to 68 galaxies for which dynam- 
ical models were created. This sample is identical to that used in 
IShettv & Cappellaril d2014lf . 


3.2 Sample Bias 

Given our strict selection criteria, our sample cannot be represen¬ 
tative of the general galaxy population at z ~ 1. To illustrate the 
effect of the selection criteria on the sample, in Fig. Q] we com¬ 
pare our galaxy samples to a larger parent sample of galaxies. This 
parent sample contains 2, 896 galaxies selected to have “secure” 
or “very secure ” redshift 0.7 < z < 0.9 in the EGS field of the 
DEEP2 survey. We compare the galaxy samples using data from 
the photometric catalogue provid ed by the DEEP2 survey, which 
is derived from iCoil et alJ 11200 41). The catalogue contains the ob¬ 
served Bab, Rab and Jab magnitudes of the galaxies calculated 
from observations made by Canada-France-Hawaii Telescope using 
the CFH12K camera. The catalogue also contains an estimate of the 
Gaussian radius of the galaxies which we use as a proxy for the size 
of the galaxies. For further details on the the catalogue, we refer the 
reader to lCoil et aD d2004h . 

In the left panel of Fig. [T] we compare the Bab — Rab colour 
of the galaxies against their absolute Rab observed band magni¬ 
tude. The plot depicts the -Rab magnitude limit of the DEEP2 sur¬ 
vey and clearly separates the blue cloud and red sequence of the 
field. The plot demonstrates that our sample have a higher Rab 
cut-off than the DEEP2 which can be explained by the S/N cut-off 
that we place on our galaxy spectra. This is further demonstrated 
in the right panel of Fig. |T] where we plot the Rab of the galaxies 
against a measure of the size of the galaxies. Due to the design of 
the survey, particularly the fact that each galaxy has a fixed expo¬ 
sure of 3 x 20 minutes and fixed aperture size, galaxies with lower 
surface brightness will be of lower quality. Hence, our S/N cut-off 
can be roughly approximated as a surface brightness cut-off. 

The left panel of Fig. Q] also demonstrates that our secondary 
sample, indicated by the red points in the plot, is a good representa¬ 
tion of the red sequence along with a few galaxies in the blue cloud. 
This is likely the consequence of the additional selection that we 
place on the secondary sample, whereby only galaxies with smooth 
light profiles and without significant star formation younger than 

1.2 Gyrs are selected. Our full sample, indicated by the red and 
the blue points, appears to represent all but the bluest galaxies over 
a certain magnitude cut-off. This under-representation may be the 
result of the visual inspection of the galaxy spectra, where the pres¬ 
ence of hydrogen emission or apparent high frequency noise make 
the spectrum fits look unreliable. 


4 METHODS 

4.1 Spectrum Quality Criteria 

We measure the mean S/N by fitting the observed galaxy spectra 
with an empirical stellar spectral library. Our choice for an empir¬ 
ical stellar library, rather than population models, to estimate the 
S/N is done to allow for maximum freedom to the optimal tem¬ 
plate combination, especially for very young ages. This avoids po¬ 
tential bias in the stellar velocity dispersion, which is also extracted 
during this step. However our results are quite insensitive to this 
choice, due to the rather low S/N of our spectra. The spectral fitting 
is done using the pPXF method and code flCappellari & Emselleml 


l2004 lt which uses a penalised maximum likelihood fitting technique 
in pixel space. Before fitting we logarithmically rebin the observed 
galaxy spectrum to 60 km s -1 per pixel. As templates we use a sub¬ 
set of empirical stellar spectra taken from the Indo-US Library of 
Coude Feed Stellar Spectra Library I Valdes et al .120041) . The subset 
of 53 templates is made such that: (i) each spectrum is gap-free and 
(ii) the subset is a good representation of the library’s atmospheric 
parameter range {T e s vs [Fe/H]). We use pPXF to fit for the veloc¬ 
ity and velocity dispersion of the galaxies using 4 additive and one 
multiplicative polynomial. During the fit of the galaxy spectra gas 
emission features are masked. The residuals of the fits are found 
to be random, suggesting that our spectral fitting is not biased by 
template mismatch. 

Using these fits, we derive the S/N of each galaxy spectra 
as the ratio of the median of the galaxy spectra and the standard 
deviations of the residuals. We remove galaxies that have S/N < 3 
from our final sample as we found the fits to their spectra unreliable. 


4.2 Recovering Star Formation History 


We derive the star formation history of our main sample of 154 
galaxi es v ia a_ full-spectrum fitting technique using the pPXF 
code dCannellari & Emselleml 120041) . For our templates, we use 
the MILES stellar population models dVazdekis et alJl20 10t. These 
modeH_ are based on the MILES empirical stellar spectra library 
dSanchez -Blazquez et al.l 2 0061) and have a resolution of 2.5 A 
FWHM i Falcon-Barroso et alJl201ll) . To ensure realistic fits of the 
galaxy spectrum with the template models, we placed certain con¬ 
straints on our template grid. We use stellar population models log¬ 
arithmically spaced into 40 age bins within the age range of 0.08- 
7.9 Gyrs to restrict the age of stellar populations to that less than 
the age of the universe at the redshift under study. We have also re¬ 
stricted our template metallicities to [M/H] of -0.4,0.0 (Solar metal- 
licit^) and 0.22. This constrain is based on results by lSchiavon et alJ 
( 2006 ) where the authors studied the absorption line strengths of 
stacked spectra of the DEEP2 sample at redshifts 0.7 < z < 1 
and concluded that these_galaxies have metallicities close to solar. 
Since ISchiavon et all d2006l) use [Fe/H] to measure their metallici¬ 
ties, we should note that for our template models the [M/H] is equiv¬ 


alent to [Fe/H] dVazdekis et al 20ld ). For this analysis we adopt 
the Salpeter IMF i SaipeteiTl955 ) for our stellar population models. 


however the derived trendsjn the SFH are insensitive to the IMF 
choice, as demonstrated in Panter et al.l |2007|) . 

In this study, we do a full-spectrum fitting for the spectral fea¬ 
tures and continuum shape of the observed galaxy spectra. This is 
done using four multiplicative polynomials to account for effects 
of reddening and calibration uncertainties in the observed galaxy 
spectra. The fractional weight assigned to each template indicates 
the fraction of the galaxy spectrum contributed by that stellar popu¬ 
lation. We do not use any additive polynomials when fitting for the 
stellar populations as these polynomials can artificially change the 
line strengths of the model templates. During the fitting with pPXF, 
the code convolves the model spectra with a Gaussian Line-Of- 
Sight Velocity Distribution (LOSVD), which is optimized to best 
fit the spectrum. 

We derive the star formation history of the galaxies by fit¬ 
ting their spectra with the above mentioned template grid while en¬ 
forcing a smoothness criterion on the distribution of the weights 
in a mass weighted manner. This smoothness criterion is enforced 
by penalising the \ 2 f° r weight distributions with non-zero sec¬ 
ond partial derivatives. In this way, for infinite regularization, the 
2-dimensional solution converges to a plane. Regularization is a 
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Test of Stellar Mass Accuracy 



Figure 2. Here we present a comparison between the stellar mass derived 
via our population modelling with the same derived bv lBundv et alj f20061) . 
The linear fitting is done via Least Trimmed Square fitting described in Cl3b 
with the code modified to fit the slope at one. The best-fitting coefficients 
for the fitted relationship, y = a + b(x — xo), are presented in the top 
left corner of the plot. The solid red line is the best fitting relation between 
the two quantities, while the dashed and dotted red lines enclose regions 
of 1A and 2.6A respectively, where A is the observed rms scatter around 
the relation. For comparison, the solid black line represents the one-to-one 
relation between the quantities. 


standa rd metho d to solve gene ral ill-posed inverse problems (e.g. 
Engl, Hanke & Neubauenll996l) like the recovery of the SFH from 
a spectrum. It is implemented in pPXF using the classic linear- 
regularization approach as in equation (19.5.10) of [Pres ] J2007h 
and is enforced via the REGUL keyword in pPXF. The value as¬ 
signed to the REGUL keyword specifies the strength of the penali¬ 
sation of the x 2 - A common guideline is to penalise the model till 
A\' 2 = X 2 - Xunregul = V% X AW, where xSnregul is the X 2 of 
the model without penalisation and TVdoF is the number of degrees 
of freedom IPres skOOl section 19.5). 

We should note that the regularization approach is not similar 
to smoothing the solution after the fit. In fact regularization does 
not prevent sharp features in the recovered star formation or metal- 
licity distribution, as long as these are required to fit the data. The 
regularization only influence the weight distribution when the so¬ 
lution is degenerate and many different and noisy ones can fit the 
data equally well. In this case the smoothest solution is preferred 
over the noisy ones. This is a general feature of all regularization 
approaches. 

For a few galaxies, we find the value of REGUL is poorly 
constrained by the Ax 2 criterion, due to the rather large noise in 
the galaxy spectra. For this reason pPXF is able to fit the spec¬ 
trum using a perfectly smooth weight distribution before Ax 2 = 
V2 x IVdoF- Hence, to be conservative, and to avoid our results 
from being driven by the regularization, we lowered the require¬ 
ment on Ax 2 from the standard one mentioned above to Ax 2 = 
x/2 x Adof/IO. 

As a consistency check, we have compared the stellar masses 
of 151_common galaxies derived in this study and that of 
iBundv et all ( 2006 ). where the authors use BRIK colour to derive 


the stellar mass for a sample of ~ 8, 000 galaxie s in th e DEEP2 
survey under the assumption of a Chabrier IMF l lChabrieill2003l) . 
We calculate our galaxy stellar mass (M+ al ) by multiplying the ab¬ 
solute luminosity of the galaxies by the (Mt/L) inferred from the 
spectral fits 


(M*/L) Sal 


E jWjL B ,j 


( 1 ) 


where Wj is the weight attached to the jth template by the reg¬ 
ularised mass-weighted fit to the galaxy spectrum, M°° gaH is the 
mass of the template that is in stars and stellar remnants, and Lb,j 
is the 5-band (Vega) luminosity of the template. Please note that all 
M/Ls derived in this study are in the B (Vega) band. We derive this 
luminosity using the MGE parametrization of the HST /ACS galaxy 
thumbnails, mentioned in Section lT3l using the galaxy redshift and 
equation (10) of Cl3b. 

The comparison between the 151 common galaxies is shown 
in Fig. [2] where we compare the two stellar masses using the 
LTSXINEFIT method Q (Cl3b) which has been modified to derive 
the best fit with a fixed slope. The solid red line in the plot repre¬ 
sents the best fitting unitary-slope relation to the points, with the 
coefficients presented in the top left corner. The A parameter is the 
observed rms scatter around the relation, while the dashed and dot¬ 
ted red lines indicate the A and 2.6A deviations from the best-fit. 
The solid black line indicates the one-to-one relation of the quan¬ 
tities. We see that a relation with slope unity is a good fit between 
the two datasets, though there is a significant offset of 0.11 dex and 
some scatter. The cause of the scatter and offset is likely the effect of 
differences in the assumed IMF betwee n the two quantities, ie this 
study assumes a Salpeter IMF while lBundv et al . (2006 ) assume a 
Chabrier IMF, and the diff erences between the ph otometry and its 
analysis. Also, the work bv lGallazzi & Bell l2009i) found that based 
on the techniques used to model the (Mt/L) of the galaxies, un¬ 
certainties and systematics of the order 0.1-0.2 dex are expected. 

Giv en that the stellar masses derived in this study and that 
of IBundv et al J ( 12006 ) are completely independent, use very dif¬ 
ferent data and approaches, and suffer from different systematics, 
we adopt A/%/2 = 0.16 dex as our very conservative errors in 
the estimation of the stellar mass of an individual galaxy via stellar 
population modelling. 

To estimate a conservative uncertainty in our derived 
(Mt/L) sal, we derive our (Mt/L) under 2 model sets and fitting 
conditions: (i) a regularised fitting to the complete model grid of 
120 single stellar population models and, (ii) a non-regularised best 
fit to the model grid limited to 40 models with solar metallicities. 
Our comparison of the two in Fig.[3]shows a linear relation between 
the two values. The observed scatter is A = 0.12 dex, which corre¬ 
sponds to an error of A/%/2 = 0.08 dex or 22% in each (Mt/L), 
assuming their errors are comparable. 

For comparison we also derived the error on the de¬ 
rived (Mt/L )sai via a more standard bootstrapping technique. 
We reshuffled the residual and repeated the determination of 
(Af*/L)sai multiple times. This yielded a median error of 0.05 dex 
in (Mt/L) s a i for the full sample, which is smaller than the error of 
0.08 dex derived above. To be conservative, and to account at least 
in part for systematic effects, in this study we adopt the intrinsic 
scatter of 0.08 dex between the two quantities as the error in the 
derived (Mt/L). 


1 Available at http://purl.org/cappellari/software 
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Robustness of Stellar (M*/L) 



log(M,/L Regul ) [M s /L 0b ] 

Figure 3. In this plot, we compare the stellar M/L derived under two dif¬ 
ferent assumptions for the stellar population of the galaxies. The compari¬ 
son between the derived quantities is done using the Least Trimmed Square 
(LTS) fitting described in Cl3b. The best-fitting coefficients for the fitted 
relationship, y = a + bx , are presented in the top left corner of the plot. The 
solid red line is the best fitting relation between the two quantities, while the 
dashed and dotted red lines enclose regions of 1A and 2.6 A respectively, 
where A is the observed rms scatter around the relation. The solid black 
line is the one-to-one relation between the quantities. The independence of 
the quantities of systematic bias indicates that our fits do not suffer from 
template mismatch. 


4.3 Photometry Parametrization 

To study the scaling relations, we measure the luminosities and size 
by parametrising the photometry of the galaxies. In this study, we 
choose the F814W photometry taken by the HST /ACS to do so. The 
need for HST photometry comes from the requirement, in the dy¬ 
namical models, of an accurate description of the tracer population 
producing the observed kinematics. This also motivates our choice 
for the F814W filter, which is as close as possible to the wavelength 
range from which the stellar kinematics is obtained. 

We parametrize the photometry of the galax- 

ies using the M ultiple-Ga ussi an Expansi on _ technique 

llEmsellem. Monnet & Bacoi] 1 19941 : ICanpellaril l2002h . This 
technique fits a series of two dimensional Gaussians, each defined 
by an amplitude, standard deviation and axial ratio, to the observed 
surface brightness of the galaxy image in a non-parametric manner, 
i.e. without the use of a predefined light profile function like 
the Sersic profiles. The advantage of this technique is that the 
convolution and deprojection (i.e. conversion of the 2D model into 
the intrinsic 3D luminosity distribution) of the MGE models is 
analytically simple. The fit is performed with the robust method 
and software (see footnote 1) by lCappellaril d2002 ). 

In this study, we use the MGE parametrisation of the photom¬ 
etry to construct our dynamical models. These models assume ax- 
isymmetry, but gala xies ca n show non-axisymmetric features, such 
as bars. IScott et al.1 d2013l Sec. 3.2.1) tried to reduce the effect of 
bars in the models by essentially ignoring them during the fits, while 
forcing the MGE to only describe the light distribution of the under¬ 
lying galaxy disk. In this study, we also adopted the same approach 
for all galaxies. 


The luminosity of the galaxies are calculated from the analytic 
sum of the Gaussian luminosities equation (10) of Cl3b, while the 
projected half-light radii ( R e ) were computed from the circularised 
Gaussians using equation (11), of C13b. For our galaxies, we also 
calculate the R ^ /Iaj , the major axis of the isophote containing half of 
the light of the galaxy as these are found to be a more robust mea¬ 
surement of the galaxy si ze due to its lower dependency on galaxy 
inclination dHopkinset alj20ld. Cl3b). We derive this value as de¬ 
scribed in section 3.3.1 of Cl3b. 

In Shetty & Cappellari d2014l) we show that, for the sub¬ 
set of 68 galaxies for which we construct dynamical models, the 
R e derived using the MGE parametrization are consistently off¬ 
set from the R e one can predict from the virial e quation A-/ V j r = 
5.0 R e tJe/G, when using the coefficient of lCappellari et all d2006l) . 
We found that those galaxies required a 0.16 dex increase in R e , for 
M v i r to match the Mjam derived via JAM models. Since we use 
the same data, at the same redshift, for all 154 galaxies, we apply 
the same fractional correction to the R e and J?;r ai of all galaxies 
of this study. Since the F814W filter approximates to the B band 
at the redshift range of our galaxies, we state all luminosity and 
photometric quantities in the B band. 


4.4 Dynamical Modelling 


To further understand the evolution of galaxies since 2 ~ 1, we 
created dynamical models for a subsample of our galaxies. In this 
study we use the Jeans Anisotropic Modelling (JAM; see footnote 
1 ) method and code of ICappellaril 1200 80, which solves the Ijeansi 
d 1922b equations under the assumptions of axisymmetry, while al¬ 
lowing for orbital anisotropy and rigorously accounting for seeing 
and aperture effects. 

The models require as input an accurate description of the 
projected light distribution of the galaxy and a parametrization 
for the unknown total mass. As the galaxy inclination is gener¬ 
ally unknown, we assume a standard inclination of 60°, the aver¬ 
age inclination for a set of random orientations. When this incli¬ 
nation is not allowed by the MGE model, we adopt the lowest al¬ 
lowed one. Here an inclination of 90° corresponds to an edge-on 
view. Another quantity required by JAM is the velocity anisotropy 
/3 Z = 1 — a^/aji, where a z and <jr are the velocity dispersion along 
the symmetry axis and along the cylindrical radius, respectively. In 
this study, we adopt (3 Z — 0.2 as this is found to be a typical value 
for ETGs in the local universe (Gerhard et al. 2001; Cappellari et al. 
120071) . Tests bv ICappellari et al] 1 20061) and LabiancheetalT 20120 
show that the derived M/L through the JAM code is weakly sen¬ 
sitive to the effects of inclination and spatially constant velocity 
anisotropy. 

The 68 galaxies in our modelling subsampl e are the same fo r 
which we created dynamical models in Shettv & Capp ellaril d2014l) . 
and we adopt the stellar velocity dispersion from that study. The 
stellar kinematics were derived from the DEEP2 spectra within a 
1 arcsec slit. In I Shettv & Cappellaril d2014b our JAM models as¬ 
sumed the total mass is distributed like the luminous one, hence 
our derived dynamical M/L incorporates the effect of stellar and 
dark matter. Here we improve our dynamical models by explicitly 
including dark matter using results from the abundance matching 
techniques. This infers a relation between the stellar and halo mass 
by matching simulated dark halo mass functions to observed galaxy 
luminosity functions, under the assumption that the most mas¬ 
sive galaxies reside in the most massive halos (e.g, KravtstwetalJ 
2004: Iconrov. Wechsler & Kravtsovll2006l: [Vale & Os trikei 20061 ; 


Vale & Qstrikeifl2008i ). 
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Here we try two different relat ions between th e ste llar and 
halo mass. The first is taken from iLeauthaud et alj d2012h . where 
the authors fit a log-normal Stellar-to-Halo mass relation (SHMR), 
with a log-normal scatter, to the results of weak galaxy-galaxy 
lensing, galaxy clustering and galaxy number density analysis of 


ixy 

the COSMOS survey l Scoville et al.ll2007h . We use equation (13) 
of ILeauthaud et alJ J20 1 2h and the parameter values for model 
SIGJAOD1 within the relevant redshift bin to derive the halo mass 
associated with a galaxy of given stellar mass. The sec ond SHMR 
we use is equation (2) of iMoster. Naab & White! ( 2013.1 . Here the 
authors use high resolution N-body simulations to derive the dark 
halo mass function at different redshifts, and match them to the ob¬ 
served galaxy stellar mass functions, hence allowing one to derive 
the halo mass for a galaxy of a given stellar mass. 

We further assume that the dark mat ter of our galaxies is dis- 
tributed spherically with an NFW profile iNavarro. Frenk & White! 
19961) . We constrain the concentration of the halo using equa¬ 
tion (12) of iKlypin, Truiillo-Gomez & Primackl d201lh . where the 
authors study the evolution of the concentration of the galaxies with 
redshift using the Bolshoi simulation. Using the above equations, 
we calculate the NFW density profile, which is fit using the one di¬ 
mensional MGE fitting code of Icappellaril d2002l- see footnote 1), 
to use in JAM. 

In practice, starting from a trial (M*/L) jam for a galaxy, we 
compute the total stellar mass by multiplication with the total lu¬ 
minosity. We then associate a corresponding dark halo mass using 
the SHMR. and a NFW profile using the halo mass versus concen¬ 
tration relation mentioned above. For given stellar and dark matter 
distributions we use JAM to predict the stellar U rms = VV 2 + er 2 , 
integrated within the adopted slit and convolved by the seeing. The 
model Vrms is matched to the observed velocity dispersion by vary¬ 
ing the only free parameter, (M*/L) of our model, using a root 
finding algorithm. In this manner we account for the dark matter 
content of the 68 galaxies in our sample, and derive their dynamical 
stellar ( M/L ) ((M*/L) jam). 


5 RESULTS 

5.1 Observed Star Formation History 

In Figj4] we have plotted the SFH of the galaxies after having col¬ 
lapsed them along the metallicity axis. The galaxies are sorted along 
the x-axis of the plot in order of increasing stellar mass as derived 
by stellar population modelling. In the top panel, we present the 
regularised SFH and we can clearly see that there is a gradual trend 
in the presented SFH with stellar mass. It appears that galaxies with 
low stellar mass tend to have a significant mass fraction of their stel¬ 
lar populations younger than 1 Gyr, unlike the typical SFH of high 
stellar mass galaxies, where the significance of this young stellar 
population weakens. 

For comparison in the bottom panel we have presented the 
unregularised mass distribution for the full-spectrum fitting of the 
galaxies. In this case the weights distribution is much noisier and 
consists of discrete peaks, as expected. But the main trend in SFH 
with galaxy mass can still be recognized. This shows that the recov¬ 
ered trend is robust. 

To further illustrate this trend, we have coadded the galaxy 
spectra into 3 bins: 10 111 < A/, Sal , 10 10 ' 75 < M* Sal < 10 11 ' 1 and 
M* al < 10 10 ' 75 . The full-spectrum fitting of these coadded spectra, 
along with their associated SFH, are shown in Fig. [5] As the S/N is 
much higher for the coadded spectrum of the galaxies, we adopt 


Trend in Star Formation History 



log(M® al /IVy-> 



Figure 4. These plots illustrate the fraction of mass represented by stellar 
populations of different ages for each galaxy. Along the y-axis are stellar 
populations of different ages, while along the x-axis are the stellar masses 
of the galaxies, sorted such that the stellar masses increases from the left to 
right. The intensity of the blue shade represents the fraction of mass of the 
galaxy composed of stellar populations of a certain age. The top panel repre¬ 
sents the SFH of the galaxies derived through the mass-weighted regularised 
full-spectrum fitting of the galaxy spectra to the entire template grid, while 
the bottom panel represents the mass distribution from a mass-weighted un¬ 
regularised full spectrum fitting to the stellar populations with only solar 
metallicities. Both plots demonstrate that galaxies with lower stellar mass 
tend to have a significant fraction of young stellar populations unlike those 
with the highest stellar masses where the bulk of their stellar populations are 
formed very early on in the history of the universe. This correlation is further 
demonstrated in Fig. [5] where we have coadded the galaxy spectra based on 
their stellar mass. The bin sizes of these stellar mass bins are represented 
here in vertical red dashed lines at stellar masses. 


the standard A^ 2 = >/2 x NboF criterion for regularization (for 
further information, please refer to Section l4.2b . 

Also, the significantly high S/N allows us to confidently fit the 
coadded spectra using the entire metallicity range of the MILES 
models, [M/H] = -1.7 to 0.22. We should point out that the [M/H] 
is n ot equivalent to [F e/H] for the low metallicity models, however 
Vazdekis et alj d201Qh state that the this doesn’t effect the galaxy 
age and metallicity estimates significantly. The same trend that we 
could see in the individual spectra appears cleaner in the stacked 
ones: (i) the most massive galaxies require a maximally old and 
nearly solar population to reproduce the observations, with little 
room for recent star formation, (ii) The galaxies with intermediate 
masses still peak at the oldest ages, but allow for a more extended 
star formation history, (iii) The lowest mass galaxies do not peak at 
the oldest ages any more and additionally show clear evidence for 
ongoing star formation. 
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SFH vs Mass trend from stacked spectra 

10 111 < AA Sal 
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Figure 5. Here we present the coadded spectra for the three stellar mass bins: 10 111 < M* al , 10 10 - 75 < M* al < 10 111 and M* al < 10 10,75 . The plots 
on the left hand side illustrate each coadded spectra (Black line) along with the best fit from pPXF (Red line) plotted on top. The green points below are the 
residuals of the fit vertically shifted to -0.5. For reference, we have also plotted the location of stellar absorption and gas emission lines within the observed 
wavelength. By eye it is clear that the quality of the coadded spectrum is far superior to that observed for each individual spectrum. In this case, the pPXF fit 
includes both stellar SSP models and gas emission line templates. For the latter, we included in the fits the Balmer emissions and the [OIIJA3727 doublet, for 
which we adopted the same kinematics. This allows us to use the full spectral range without the need for masking the emission lines. Along the right hand side, 
we have plotted the mass fraction of the fitted stellar population models of various ages and metallicities. We clearly see a trend between the star formation 
history of the galaxies along with the stellar mass of the galaxies. The coadded spectra of the highest masses indicate that these galaxies are likely to form the 
bulk of their stars in the early epochs of the universe and undergo a passive evolution thereafter. In the case of the intermediate masses, the galaxies appear to 
have an extended SFH while the lowest mass spectra require recent star formation, along with a very old population, to match the observed spectra. 
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The well known age-metallicity degeneracy of stellar popula¬ 
tion models implies that the effect on a galaxy spectrum due to a 
decrea se in its mean age can be negated by an increase of its metal- 
licity (Worthevlll994l : lRenzinil l2006h . This degeneracy is reduced 
by the use of full-spectrum fi tting , which tries to reproduce a large 
number of spectral features (Conrovlf2013l ). but the general trend is 
still present in our recovered age and metallicity distributions. An 
example of this can be seen in the bottom right panel of Fig. [5] 
The ridge of equally-large weights going from an age of 2 Gyr with 
nearly solar metallicity, to an age of 7 Gyr with the lowest metal¬ 
licity is likely due to the degeneracy, combined to our relatively 
narrow wavelength range. In Table. 1, we present the relevant stel¬ 
lar modelling quantities of the full galaxy sample of the study. The 
complete version of this table is available online with the published 
article. 

In this study, we calculate the mean age of the galaxies from 
the SFH as follows: 

/, / * ,, J2 j Wjlog(Age j ) 

(log(Age)) = „ -, (2) 

2^j Wj 

where Wj is the weight associated to the j th template during the 
regularised luminosity-weighted full-spectrum fit and Age^ is the 
corresponding age. To obtain luminosity-weighted quantities one 
can simply normalize all the SSPs to the same mean flux within the 
fitted wavelength range, before the pPXF fit. Of course the lumi¬ 
nosity weighting refer to the fitted wavelength, which in our case 
lies between the U and B band. When the SSP spectra used in the 
pPXF fit have the flux as given by the models, namely correspond¬ 
ing to a unitary stellar mass, this equation provides mass-weighted 
quantities. 

In Fig. [6] we present the luminosity-weighted ages of the 
galaxies against their respective stellar masses. The black points 
represent the full sample of 154 galaxies of this study, while the red 
triangles represent the luminosity-weighted ages of the entire sam¬ 
ple of ~ 70 galaxies of G14, at similar redshift. On the bottom right 
of the figure, we plot a representative error bar for the 154 galaxies 
in this study. The error on the stellar mass is based on the scatter 
in Fig. [2 ie 0.16 dex. To estimate the error on the ages we boot¬ 
strap the galaxy spectra to derive their age multiple times. From 
this distribution of ages for each galaxy, we calculate robust stan¬ 
dard deviations and take the median of these as the representative 
error. This bootstrapping is done using only the templates with solar 
metallicities and without the regularisation of the template weights. 

The two sample appear to largely overlap, however our sample 
includes a larger number of low-mass and correspondingly younger 
galaxies. The solid black line is the resulting best-fit to our galaxy 
sample using a similar formula as G14: 


Age b (M b y 

Age(AT*) \Mj ’ 


(3) 


which describes an age trend as a function of stellar mass M* fol¬ 
lowing a power-law Age oc M* at small masses and reaching the 
asymptotically constant value Age b above the characteristic mass 
Mb, with a smooth transition in between. The best-fitting coeffi¬ 
cients for the luminosity-weighted and mass-weighted (not shown) 
ages of our 154 galaxies are presented in Table. 2. 

In Fig. [6] we also present two local samples. The median 
luminosity-weighted age trend of the sample of both spirals and 
ETGs from 'Gallazzi et alj (2005) is represented by the red solid 
line, with the red band enclosing the 16 th and 84 th percentiles of 
the age distribution. The blue contour lines represents the probabil¬ 
ity density distribution of the single stellar population ages (similar 


Stellar Mass vs. Age 



Figure 6. In this plot, we present the observed trend between the regu¬ 
larised luminosity-weighted stellar ages of galaxies and their stellar masses. 
The black points represent our sample of 154 galaxies at redshift 0.7 < 
£ < 0.9 with the solid black line representing the best fit to the trend us¬ 
ing Eq. [3] A representative error bar for the galaxies is presented in the 
bottom right of the plot. The red tria ngles in the plot are the luminosity- 
weighted ages of lGallazzi et alJ (12014b . The plot demonstrates the trend be¬ 
tween our luminosity-weighted ages and the stellar mass of the galaxies. In 
this plot we also present two samples of the local universe for reference. 
The solid red line represents the best fitting relations for the galaxies of 
I Gallazzi et alJ J2006h with the red band enclosing the 16 th and 84 th per¬ 
centiles of the distribution. The blue contours represent the probability den¬ 
sity of the ATLAS 3D sample. The dash-dotted red line represents the trend 
expected from the local galaxies of lGallazzi et alj d2006h at z ~ 0.8 if these 
galaxy have evolved passively since then. 


to luminosity-weighted ag es) of the vo lume-limited 260 ETGs of 
the ATLAS 3D s urvey JCappellari et alj 1201 lah , as determined by 
iMcDermid et al.1 (20151 Hereafter MD15). This was calculated us¬ 
ing the Epanechnikov kernel, with a smoothing template that was 
optimised for each dimension, based on lSilvermanl (19861) . The dot- 
dashed red line in the plot represents for reference the passive back 
evolution of the trend observed in lGallazzi et alj (20051) . These two 
samples provide a local universe reference for our results, and given 
the differences in the two samples selection give a sense of the pos¬ 
sible ranges in the local trends. Our most massive galaxies are not 
consistent with passively e volvi ng into the corresponding sample of 
spirals and ETGs of the lGallazzi et alJ (2005 ) sample, but they are 
consistent with evolving passively into the massive ETGs of MD15. 
This is not only true according to their mean ages but, most impor¬ 
tantly, it is also true also for the actual SFH, which in both cases 
shows the galaxies essentially consists of a maximally old stellar 
population. 

The results presented in Fig[4] and Fig[5] are in line with the 
results of G14, where the authors predict that the SFH of galaxies 
has to be a function of their stellar mass. We note that since our 
analysis is based on the spectra of galaxies taken within the central 
1” of the galaxies, these results are relevant only to that region. Our 
analysis is unable to give any direct insight into the evolution of the 
outer regions of these galaxies. 

By directly studying the SFH of galaxies at half the age of 
the universe, our study serves as a stringent test to the scenario of 
galaxy formation and evolution that was presented in the Introduc- 
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Table 1. Results of Stellar Population Modelling. The table below is a partial representation of the complete set, which is 
available as part of the online material of the paper on MNRAS. 


DEEP2 ID 

Redshift 

Re 

R f aj 

M b 

(M,/L)g al 

Mass- 

Luminosity- 

log(M Sa!) 






(Regul.) 

weighted Age 

weighted Age 




(kpc) 

(kpc) 

(Mag.) 

(Mq/Lq) 

(Gyrs) 

(Gyrs) 

(M©) 

U) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

11044628 

0.824 

9.7 

10.0 

-22.06 

0.5 

1.65 

0.28 

10.686 

11050845 

0.840 

4.4 

5.2 

-21.76 

3.5 

4.85 

4.57 

11.429 

11051246 

0.777 

4.4 

4.6 

-21.18 

0.9 

2.47 

0.71 

10.592 

12003915 

0.749 

3.2 

5.1 

-20.89 

0.7 

2.52 

0.34 

10.391 

12004067 

0.748 

9.4 

10.2 

-21.76 

0.8 

2.26 

0.54 

10.784 


Column (1) : DEEP2 galaxy identifier. Column (2) : DEEP2 estimated redshift. Column (3) : Effective radii of the galaxies 
derived analytically from the circularised MGE model. This R e has been corrected for the underestimation that was shown 
in lshettv & Cappellaril 12014 3. Column (4) : Major axis of the isophote containing half of the observed light of the galaxy. 
Similar to the R e , this value has also been corrected for the underestimation. Column (5) : Absolute B (Vega) band magnitude 
of the galaxies derived analytically through MGE models of the F814W photometry. Column (6): The mass-weighted M*/Lb 
derived through the regularised full-spectrum fit using the entire template grid. In this study, we adopt the intrinsic scatter 
between the (M*/L) derived through the different methods and template grids as their error, i.e. 0.08 dex. Column (7) : 
Mass-weighted age of the galaxies calculated from the regularised full-spectra tit using the entire template grid. Column (8) : 
Luminosity-weighted age of the galaxies calculated from a regularised full-spectra fit using the entire template grid. Column 
(9): The log of the stellar mass, in unit s of A/,-, , derived using the values in Column (7) and Column (5). Adopting the intrinsic 
scatter between our value and that of lBundv et alJ (2006), 0.16 dex is a good approximation to the error in this value. 


Table 2. Best fitting parameters for equation]?] 


Technique 

log Age 6 

log M b / Mq 

7 

Mass-weighted Ages 

9.70 

10.66 

1.15 

Luminosity-weighted Ages 

9.83 

11.35 

1.18 


tion of this paper. Work by iOnodera et alJ J2012I) . where the au¬ 
thors use stacked spectrum of 17 passively evolving galaxies at 
redshifts of 2 > 1.4, demonstrated that those massive galaxies 
ha d been quenched sin ce 2 ~ 2.5. This is consistent wi t h wor k 
by ljuneau et alJ l l2005l) and iBehroozi. Wechsler & Conrovl d2013l) . 
where the authors use observational data covering a large redshift 
range, to show that the star formation rate peaks at redshift of ~ 3 
for the most massive galaxies and that this peak evolves to lower 
redshifts for lower mass galaxies. 

This trend in the star formation rate of the galaxies has been 
known for about a decade from observations of the fossil record 
in the local Unive rse (e.g. lHeavens et all2004l : lThomas et alj|2005l : 
IPanter et~aL 2007, MD15). However the local studies require a sig¬ 
nificant extrapolation to describe the SFH over the full evolution 
of the Universe. While our study directly recovers the SFH at an 
epoch when the Universe was half its current age, it is reassuring to 
find that the extrapolated local SFH trends are consistent with what 
we observe. Our result further strengthens these results and the pic¬ 
ture of galaxy evolution that they support, by demonstrating that the 
bulk of stars of the massive galaxies had formed in the early epochs 
of the Universe and have been evolving passively, while low mass 
galaxies were/are still forming significant mass fractions of their 
stellar populations till 2 ~ 0.8. 


5.2 Mass-Size Distribution 

In Cl3b, the authors demonstrate that by replacing the luminosity 
with mass, and the circularised R e with 7Z™ aj , the Fundamental 


Plane evolves into a Mass Plane, which follows the virial equation. 
This shows that the Fundamental Plane is due to virial equilibrium, 
combined with a smooth M/L variation with the galaxy velocity 
dispersion a. Most of the information on galaxy evolution is con¬ 
tained in the non edge-on views^of the p lane , e.g. the Mass-size 
and Mass-cr distribution. Cappe llari et all 12013a. Hereafter C13a) 
shows that galaxy properties, namely the stellar age, colour, M/L, 
IMF and gas content, all best follow lines of constant a, which 
they demonstrate traces the bulge mass fraction. MD15 additionally 
showed that metallicity and non-solar abundance are also driven by 
a. Below we present the trend of the stellar population in the mass 
plane of galaxies at 2 ~ 1 and compare it with that observed in the 
local universe. 

Using the physical parameters derived in this study we are able 
to study the mass plane of our galaxies which are presented in Fig[7] 
In the left panel of Fig. [7] we plot the stellar mass, derived through 
stellar population modelling, and the galaxy size for the 154 galax¬ 
ies in our full sample. In the lower right corner, we plot a represen¬ 
tative error bar for the derived quantities. The error on the stellar 
mass derived using stellar populations modelling is based on the 
scatter observed in Fig. [2] ie 0.16 dex. The error on the galaxy s ize 
is based on our comparisons with lFemandez Lorenzo et al] 1201 ll) . 
where the authors fit a de Vaucoleurs stellar profile to I-band pho¬ 
tometry of 40 galaxies present in our sample. We compare the ef¬ 
fective radii calculated by these authors with that calculated using 
the MGE models of this study and derive an error of 0.05 dex on 
the galaxy size. The colour of the symbols in the two plots repre¬ 
sents the luminosity-weighted age (Gyrs) of the_stellar population 
in the galaxy. This was LOESS smoothed Icievelandil979 ) using 
implementation (see footnote 1) of Cl3b. This plot can be directly 
compared to fig. 1 and fig. 3 of C13a and to fig. 7 of MD15, where 
the authors respectively show that for the ATLAS 30 volume lim¬ 
ited sample of ETGs the observed stellar population indicators and 
stellar age of the galaxies correlate with their velocity dispersion. 
Our plot suggests that this trend is in place by 2 ~ 1. 

The solid red line in the left panel of Fig. [7] is the Zone of 
Exclusion (ZOE) for the nearby galaxies as prescribed by C13a, 
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Luminosity-weighted Age (Gyr) Luminosity-weighted Age (Gyr) 



Figure 7. Here we present in the left panel the mass size plot for the 154 galaxies in our full galaxy sample. The mass in the plot is derived from the mass 
weighted regularised stellar population modelling while the galaxy size is represented by ] . The solid red line of the plot defines the Zone of Exclusion 
(ZOE) determined by the ATLAS 3D group (C13a), beyond which the probability distribution of galaxies is negligible. The dashed blue lines are lines of constant 
velocity dispersion. The colour trend observed in the plot is the LOESS smoothed (Cl3b) luminosity weighted age. This plot looks similar to its dynamical 
variant in the local universe (See fig. 6 and fig. 7 of MD15), where the authors find that the age of the stellar population within a galaxy is correlated to the 
its velocity dispersion. The left panel also demonstrates that some quiescent galaxies are on the other side of the ZOE and would require to grow by a factor 
of ~ 1.5 to be consistent with local observations. In the right hand panel, we present the dynamical stellar mass, derived using the (M*/L)jam. against the 
velocity dispersion for our secondary galaxy sample, which is equivalent to the Faber-Jackson relationship. In this plot, the dashed blue lines depict lines of 
constant size. Since the quantities of this plot is derived in a manner independent of those in the left panel, this plot illustrates the robustness of the results 
discussed. In the lower right comer of both panels, we present a representative error bar for the quantities. 


beyond which the probability of finding galaxies within a volume- 
limited sample is ~ 1%. The left panel demonstrates that some 
galaxies in the full sample would lie a factor of ~ 1.5 below 
the ZOE in the local universe. The robustness of this result is 
demonstrated in the right panel of Fig. |7] where we have plot¬ 
ted the dynamical stellar mass against the stellar velocity disper¬ 
sion for the secondary galaxy sample. This plot is equivalent to 
the Faber-Jackson relationship with the galaxy luminosity substi¬ 
tuted for the dynamical stellar mass derived by the multiplication of 
(M* / L ) jam with the total luminosity of the galaxies. The solid red 
line in this plot is the ZOE translated into this Faber-Jackson plot via 
the Virial equation. Despite the stark contrast between the method¬ 
ology used to derive the physical parameters in the two plots, both 
plots present a consistent picture. Considering that the ATLAS 3D 
sample is volume-limited and thus provides a representative refer¬ 
ence for the local population, while our 2 ~ 1 sample is biased 
towards the most dense objects (Sec. 13.2^ . this implies that the fac¬ 
tor of 1.5 provides a stringent upper limit to the to the growth of 
galaxies since 2 ~ 1. To aid the reader in accurately interpreting 
this plot, a representative error is presented. This error bar is based 
on the median error of the velocity dispersion, derive d via boot¬ 
strapping (see section 3.1 of lShettv & Cappellaril (2014f)), and the 
propagation of this error onto the derived dynamical mass via eq. 8 
of lCappellaril d2008l) . 

Our upper limit to the size growth of the most compact 
g alaxies since 2 ~ 1 is consistent with previous works such as 
INewman et al.l d2012h . where the authors fit the observed mass-size 


relation for galaxies in bins of redshift. Their fits suggest that a 
galaxy of mass M* = 10 11 Mq within their redshift bin of 0.4 
to 1 grow by a factor of 1.2 by redshift of 0.06. In the work pre¬ 
sented by Ivan der Wei et alJ (2014), the authors use 3D-HST and 
CANDELS data to study the evolution of the mass-size distribution 
of 30,958 galaxies with stellar masses > It) 9 Mq and 0 < 2 < 3. 
Their comparison of the size evolution of the median sizes of galax¬ 
ies from 2 = 0.75 to 2 ~ 0 demonstrates a size growth of 
~ 1.4, consistent with our findings. There have been recent sug¬ 
gestions that part of the observed size growth in quiescent galax¬ 
ies may be due a progenitor bias, i.e. the increase in the average 
size of the quiescent galaxy population with cosmic time is due 
to the inclusion of larger recentl y quenched galaxie s into the qui- 


escen t galaxy population (e.g. Poggi anti et alj|2013 : Carolloet al 
2013) . However, the results of re cent work (e.g. iKrogager et al 


2014: Belli. Newman & Ellisi201S) have shown that the progenitor 


bias isn’t sufficient to accommodate the size growth that has been 
observed since 2 ~ 2.5. 


5.3 Dark Matter Fraction and IMF 

In Ishettv & Cappellaril d2014h . we had modelled the stellar kine¬ 
matics of the galaxies assuming that mass-follows-light and us¬ 
ing the total M/L to account for the dark matter contribution. 
When making our inferences on the IMF of those galaxies we 
made the empirically-motivated assumption that the dark mat¬ 
ter fractions within their central regions were negligible. The 
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Dark Matter Fractions 



Figure 8. Here we present the dark matter fraction within 1 R e for the 68 
galaxies for whom we have created dynamical models. The dark matter 
fraction is the dark matter enclosed within a radius of 1 Re divided by the 
total mass (i.e. stellar + dark matter) enclosed by that radius. In the top 
panel, the dark matter haloes are assigned based on the SMHR derived by 
iLeauthaud et alj 42012b while the SMHR by Moster et al (2013) is used to 
derive the dark matter haloes for the bottom panel. Despite the difference 
in the absolute value of the dark matter fraction, both results demonstrate a 
low dark matter fraction for the inner regions of massive galaxies. 


basis for our assump ti on w ere mainly the results of Cl3b and 
iHilz. Naab & Ostrikeil d2013l) . In Cl3b, the authors were able to 
derive the dark matter fraction of their 260 local ETGs using 
integral-field data and dynamical modelling, and accounting for the 
dark matter via multiple dark matter profiles. Their results robustly 
proved that the dark matter fraction of the central regions of mas- 
si ve ETGs is <30 per cent. Furt hermore, galaxy merger simulations 
of lHilz. Naab & Ostrikeil d2013h find that the dark matter fraction of 
galaxies would increase with time, i.e. with decreasing redshift. 

A robust test of the assumption would be a direct measure¬ 
ment of the dark matter fractions in these galaxies, however this is 
not feasible from the stellar dynamics without spatially-resolved in¬ 
formation, which is still difficult to obtain with current instruments. 
Therefore, to test this assumption, we place our galaxies within dark 
matter haloes predicted by abundance matching techniques. Fur¬ 
thermore, to get robust results on the dark matter fractions of the 
galaxies, we have derived the dark matter fractions for our galaxies 
using two independent SHMRs (Sec. l4.4l l. 

The results of the dynamical modelling are presented in Ta¬ 
ble. 3 and illustrated by Fig. [8] The plots demonstrate that, to be 
consistent with the observed surface brightness and velocity disper¬ 
sion, as well as with the halo matching results, massive galaxies at 
this redshift must contain low dark matter fraction fr>M{Re) within 
a sphere of radius 1 R e . The inferred median value for our sample 


log a (km/s) 

1.80 2.14 2.47 


■ ■ I 



log(M*/L) JAM [M 0 /L 0b ] 


Figure 9. Here we present the IMF normalisation for the 68 galax¬ 
ies \yith dyna mical models. This is an updated version of the fig. 4 
in IShettv & Cappellaril 42014b the difference being that here we use the 
(M*/L)jam derived from dynamical models, instead of the total one. 
Along the x-axis, we have plotted the (M*/L) derived via dynamical mod¬ 
elling while the y-axis represents the stellar (M* / L ) derived via stellar pop¬ 
ulation modelling, assuming a Salpeter IMF. The colour trend shown here 
is the LOESS smoothed (Cl3b) central velocity dispersion of the galaxies. 
Representative errors on the (M*/ L) is presen t ed in the top left comer of 
the plot, based on Fig. [2] and eq. 8 o fCappellarl 12008 4 


of 68 galaxies is fr>M(Re) = 9 and 13 per cent, while the 90th 
percentile is /dm (-Re) = 21 and 34 per cent, using the SHMR of 
ILeauthaud et alj d2012l) and lMoster, Naab & White d2013l) respec¬ 
tively. 

In the two panels of Fig. [8] the green line represents the ro¬ 
bust parabolic fit to the dark matter fractions of the 260 galaxies 
(Cl3b) derived using SHMR relation of lMoster et all bold) . while 
additionally requiring the models to best fit the integral-field kine¬ 
matics and photomet ry. This is si milar to our dark matter fraction 
results derived using lMoster. Naab & Whitel 120131) . Comparing our 
results to that of ATLAS 30 , the dark matter fraction appears to have 
increased in the last ~ 7 Gyrs. The lower dark matter fraction de¬ 
rived at higher redshif ts is consistent with th e numerical simulation 
prediction of IHilz. N aab & Ostrikeil 1 20 3l ) and is also consistent 
with recent results bv lBeifiori et alj a2014f) . which found that the ra¬ 
tio of the dynamical mass (dark matter + stellar mass) to the stellar 
mass of the galaxies reduces with increasing redshift. However, we 
do note that this trend may be in part influenced by the differences 
in the sample selection. 

In the work presented in lShettv & Cappeliaril ( 2014b . the IMF 
normalisation of the 68 dynamically modelled galaxies was shown 
in fig. 4. The plot illustrated that the IMF normalisation for mas¬ 
sive ETGs at 2 ~ 1 was Salpeter-like, however a caveat of the 
work was the assumption that the dark matter fractions in the central 
regions of high redshift galaxies was not significant. In this study 
we include dark matter explicitly and in Fig.[9]present an updated 
version of the plot, with the dark matter contribution explicitly re¬ 
moved. This confirms the assumption that dark matter fraction has 
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Table 3. Results of Dynamical Modelling. The table below is a partial representation of 
the complete set, which is available as part of the online material of the paper on MNRAS. 


DEEP2 ID 

(1) 

a 

(km s -1 ) 
(2) 

6(7 

(km s 1 ) 

(3) 

{M*/L)j AM 

(Mq/Lq) 

(4) 

log(M } AM ) 
(M®) 

(5) 

fDM 

(6) 

11050845 

298.2 

28.7 

5.74 

11.644 

0.088 

12004136 

202.4 

21.7 

2.80 

11.275 

0.178 

12004516 

176.4 

28.6 

2.56 

11.164 

0.146 

12008254 

250.7 

33.3 

4.56 

11.444 

0.103 

12008360 

185.2 

31.1 

2.54 

10.974 

0.035 

12008441 

257.7 

64.3 

8.53 

11.646 

0.141 

12011900 

318.5 

57.1 

9.30 

11.499 

0.042 


Column (1) : DEEP2 galaxy identifier. Column (2) : Velocity dispersion as measured 
by full-spectrum fitting of the galaxy spectra. Column (3) : Error on the derived velocity 
dispersion of the galaxy through a bootstrapping technique. Column (4): B-band (M* / L) 
derived th rough dynamical modelling, where the dark matter halo was assigned using the 
SHMR of iLeauthaud et all d2012t) . Column (5) : The log of the stellar mass, in units of 
Mg, of the galaxy derived using the absolute luminosity (Column (5) of Table. 1) and the 
M/L of Column (4). Column (6) : fraction of' dark matter within 1 R e of the galaxy. This 
was derived by using the SHMR of lLeauthaud et al.U2012l) . 


a negligible effect on the IMF normalisation of the galaxies, i.e. the 
average IMF normalisation of massive ETGs at redshift of ~ 1 is 
still consistent with a Salpeter IMF. This result is independent of 
the SHMR used. 

This confirmation of our lshettv & Cappellari) 1 2014|l result is 
also in agreement with recent work bv Isoimenfeld et al.1 J2015h . 
where the authors have used a joint lensing and stellar dynamics 
modelling of massive galaxies out t oz« 0.8, including about ten 
galaxies within our resdhift range. The authors found that the stellar 
IMF nor malisation was close to Salpet er IMF for M* = 11.5. Also, 
work by iMartfn-Navarro etaL 1 2015 :1. using IMF sensitive spec¬ 
tral features, determined that more massive galaxies have a bottom 
heavier IMF at redshifts between 0.9 and 1.5. The Salpeter normal¬ 
ization of the IMF is also consistent with the non-universality of the 
IMF r eported in the nearby U niverse (e.g. Ivan Dokkum & Conrovl 
l2010l: Cappella ri et alj |2012|) . In fact, if the centres of massive 
galaxies at a ~ 1 passively evolve into the massive ETGs popu¬ 
lation, their IMF is expected to have a Salpeter normalization as 
observed locally j Auger et alJl201(jl : IConrov & van PokkuidEoTl 


observ ed locally (Auger 

ICappellari et alJ2013ai .~ 


6 SUMMARY 

In this study, we have determined the SFH of the galaxies at redshift 
0.7 < z < 0.9 in a non-parametric manner, using full-spectrum fit¬ 
ting. From a parent sample of 2, 896 galaxies from the DEEP2 sur¬ 
vey, we apply strict quality selection criteria to extract 154 galax¬ 
ies with 10 10 < M» < 10 12 Mq for which we derive our SFHs. 
For a subsample of 68 galaxies, with M* > 10 11 Mq, we addi¬ 
tionally construct dynamical models. Due to these selection crite¬ 
ria, our galaxy samples have a higher magnitude cut-off than the 
DEEP2 survey and preferentially selects galaxies with high surface 
brightness. However, we show that the full galaxy sample is repre¬ 
sentative of the red sequence and all but the bluest galaxies of the 
blue cloud, while the secondary sample mainly consists of galaxies 
from the red sequence. 

The derived SFH for the full galaxy sample indicates that the 
most massive galaxies formed the bulk of their stars in the early 


epoch of the universe unlike low mass galaxies which are forming 
stars at a significant rate even at redshift of ~ 1. This is qualitatively 
consistent with previous fossil record studies in the local universe, 
where authors have found that the star formation rate of the most 
massive galaxies peaked early in the age of the universe, and hence 
provides a robust and independent test of these results and the nar¬ 
rative of the formation and evolution of galaxies that these results 
have produced. Our results based on Fig.[4]and Fig.[5]demonstrate 
that the difference between the SFH of galaxies evolves gradually 
and is a function of stellar mass. 

We study the distribution of galaxy ages on the mass-size dia¬ 
gram. This demonstrates that the velocity-dispersion dependence in 
the age of the stellar populations of the central regions of the galax¬ 
ies was already in place by 2 ~ 1. We also place an upper limit 
of a factor of ~ 1.5 to the size growth of individual galaxies since 
z ~ 1, in agreement with other studies. 

Finally, using the dynamical models of the 68 galaxies in our 
secondary sample, which account for the dark matter in the galax¬ 
ies using results from abundance matching techniques, we measure 
a median dark matter fraction /dm (-Re) ~ 10 per cent per cent, 
within a sphere of radius R e , for the most massive galaxies, with 
small variations depending on the adopted Stellar-Halo mass rela¬ 
tion. Comparing this to the dark matter fraction determined locally, 
we find that the dark matter fraction of galaxies has marginally in¬ 
creased in the last 8 Gyrs, but is otherwise insig nificant. This result 
confirm the study of lShettv & Cappellari II2014 9 stating that the av¬ 
erage IMF normalisation of the most massive galaxies is on average 
consistent with a Salpeter IMF. 
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